Engineering mixed states in a degenerate four-state system 
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A method is proposed for preparing any pure and a wide class of mixed quantum states in the 
decoherence-free ground-state subspace of a degenerate multilevel lambda system. The scheme is a 
combination of optical pumping and a series of coherent excitation processes, and for a given pulse 
sequence the same final state is obtained regardless of the initial state of the system. The method 
is robust with respect to the fluctuation of the pulse areas, like in adiabatic methods, however, the 
field amplitude can be adjusted in a larger range. 
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Controlling the quantum state of degenerate quantum 
systems have drown much attention recently. This field 
has developed independently in several different series of 
studies: Among the numerous adiabatic passage tech- 
niques one of the most well-known method is the 
stimulated Raman adiabatic passage (STIRAP) |2(. The 
STIRAP can be used not only for transferring population 
between two quantum states using crafted laser pulses, 
but it has been utilized to create coherent superpositions 
in three- and four-level systems 0, E113 , to prepare max- 
imally coherent superposition states ja| and arbitrary co- 
herent superpositions in A-state degenerate sys- 
tems. The applicability of the STIRAP method is limited 
by constraints on the field amplitudes @>E1- 

The other field that developed toward the quantum 
control of degenerate systems is termed "coherent con- 
trol" that uses several interfering pathways in the quan- 
tum system to transfer selectively population from an 
initial state to a target one ^lj. Merging this tech- 
nique with the STIRAP method led to the mapping 
of wave-packets between vibrational potential surfaces 
in molecules for the non-degenerate |l2j and degenerate 
[HQ cases. 

The above mentioned control processes have great im- 
portance in many areas of quantum-information process- 
ing (QIP), involving quantum computing, cryptography 
and teleportation |lj| . In general, mixed states cannot be 
created with coherent state preparation methods. How- 
ever, for several QIP problems it is essential to develop 
quantum state preparation techniques which are capable 
to prepare not only pure, but mixed states of the system 
as well 

In optical pumping processes , the final state of the 
system is largely independent of its initial state, how- 
ever, the efficiency is small 20]. On the other hand, in 
the coherent state-preparation methods the final state 
depends on the initial state of the system, but the effi- 
ciency can be nearly unity 0,0] • I n this letter we consider 
a novel concept for quantum-state preparation, which is 
a combination of optical pumping and coherent excita- 



tion processes, exhibiting only the advantageous prop- 
erties of the two schemes, and capable to prepare not 
only pure but prescribed mixed states of the system too. 
The unique features of our method compared to other 
state-preparation methods are the following: it is simul- 
taneously (i) robust, (ii) the final state is independent of 
the initial state of the system, (iii) the state is prepared 
in a decoherence-free subspace, (iv) the choice of the ex- 
citation field amplitude is quite arbitrary. The method 
can be implemented in multilevel lambda systems. For 
concreteness, let us consider the four-state system shown 
in Fig. ^ there are three degenerate ground states and 
a single excited state coupled by an elliptically polarized 
coherent laser pulse. The ground states \g q ) (q = — , n, +) 
are assumed to be the magnetic sublevels of a J g = 1 an- 
gular momentum state, whereas the excited state |e) has 
J e = 0. The three polarization components of the cou- 
pling field, denoted by E q with (q = — , 7r, +), share the 
same time-dependence, but they can have different peak 
amplitudes and phases 
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where the parameters 6, cp describe the polarization of 
the pulses, £ is the absolute phase of the pulse, and the 
phases of the <r+ and <r_ components relative to the it 
component are /i+ and respectively. The excited 
state |e) decays with spontaneous emission to the ground 
states with a rate of 7i n , and it may decay to states other 
than the three ground states with a rate of 7 cx t as well. 
When decay out of the ground-state space occurs, a re- 
pumping process with a rate of R p is switched on in order 
to compensate for the population loss. We show that by 
using a predetermined sequence of pulses we can create 
any prescribed pure or a wide class of mixed final states 
in the ground-state space, starting from any initial state. 

The Master equation describing the time evolution of 
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FIG. 1: (Color Online) The coupling scheme for our state- 
engineering procedure: the lower states are the magnetic sub- 
levels of a J — 1 angular momentum state which are coupled 
by <j± and n polarized pulses to a single excited state. The 
excited state decays with a rate ji n into the lower states, and 
it may decay out of the system with a rate 7 ex t ■ For non-zero 
7cxt a repumping is switched on with a rate R p . 



the system is given by 



d - 1 r ~ 

dt g = Th [H ' e] 
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where the Hamiltonian H reads 



H = -(Q_\g_)( e \ + QM(e\ 



HA\e)(e\, 



^+l.9+)<e|+h.c.) 

(3) 



where the Rabi frequencies are = 
ifie 4 «e v - sin0 sin 95, Q„ = -±fte ie cos0, Q+ = 
|fie^e v+ sin0cos<p, with ftfi = d ge £ . The step opera- 
tors are defined as L q = -±=\g q )(e\ and L e = |e)(e|. The 
other symbols are defined above. 

The Hamiltonian of Eq. © has two uncoupled eigen- 

states |4>^) [2l|, i.e. they are decoupled from the ex- 
ternal driving field, H\$<$) = for I = 1,2. They read 
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\g q ), 1 = 1,2, 
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where the unit vectors «A are given by 
= [e 1 ^- cos 9 sxsxtp, sin 6*, e v + cos 9 cos <p] T and 



n 



(2) 



cos</?, 0, e _v_ sin<p] T , and the field 
parameters fi±, ip, 9 are defined in Eq. Q). These states 
are dark states, because they do not have a component 
in the excited state |22| . The Hamiltonian has two other 
eigenstates with non-zero eigenvalues, they are called 
bright states, because they do have a component in the 
excited state [22^ . 

Let us assume that we have some initial state g- m de- 
fined in the decoherence-free, ground-state space. The 
Hamiltonian part of the Master equation J5J drives the 
bright components of this state to the excited state back 
and forth via Rabi oscillations. As the excited state gets 



populated the spontaneous emission will interrupt the 
Hamiltonian dynamics and the system falls back into the 
ground-state space or some other external states become 
populated. As a result, the two dark states become more 
and more populated, even though they are decoupled 
from the external driving field. When there is no de- 
cay into external states or the decay out of the four-state 
system is compensated by an incoherent repumping pro- 
cess, the system will relax into the dark subspace of the 
Hamiltonian. Then, the state of the system is given by 



2b rt =P w l*S , ><*g , l+P (a) l*g , X*g , l 



(5) 



where the coefficients p^ depend on the applied pulse 
and the initial state as well. It is important to note that 
this output state is independent of the pulse amplitude 
£, it depends only on the polarization and relative phases 
of the three components of the field of Eq. (JIJ . 

In the ground-state space the state l^^^) orthog- 
onal to the dark states of Eq. J3J is = 
e 1 ^- sin#sin</? \g~)— cos 9 |(7 7r )+e v + sin 9 cos <p \g+) . This 
vector can point to anywhere in the three-dimensional 
ground-state space, depending on the laser-field parame- 
ters. Consequently, it is possible to choose the laser-field 
parameters, so that any two linearly independent state 
vectors and \ip2) of the three-dimensional ground- 
state space lay in the dark subspace. Therefore, in prin- 
ciple there exists a pulse-sequence, such that any pre- 
scribed final state of the form 



Qf = Pi\ipi){ipi 



■pf ] H2)m 



(6) 



can be obtained. The state Qf can be either a pure state 
if one of the coefficients pp vanishes, or a mixed state if 
both of them are non-zero. 

We have two tasks now: (i) to find how an initial state 
gin transforms when the pulses Eq. Q are adjusted to 
a certain value; (ii) to find the pulse-sequence that steer 
the state of the system to a prescribed final state defined 
by Eq. ©. 

For convenience, the linear space of the density op- 
erators {g} is represented by vectors {r} with compo- 
nents r4(j_n + j = {g)i,j, where {g)ij is the matrix el- 
ement of the density operator g in the ordered basis 
{!#-), I Sir), \g+), |e)} ■ The scalar product of vectors is de- 



fined as (rW|r( 2 )) = £ s 
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Master equation (|3J in this representation takes the form 
4zT = Mr + d, where the matrix M describes the linear 
part of the Master equation J5J), and d corresponds to 
the constant term R p L e in the incoherent repumping of 
the excited state. In this letter we are going to consider 
two cases: 

(a) The case "f ex t = R p = : In this case the Master- 
equation is homogeneous in r, and d is zero. The relax- 
ation of the system into its final state can be described by 
those left- and right-hand eigenvectors (denoted by 
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and respectively) of the matrix M, which belong to 
the eigenvalue zero 



Mr^ } = , r 



(fc)T M = 0, (7) 



and are orthonormal (r^ \r^) = Ski- The left- and 
right-hand zero subspaces of M are four-dimensional and 
they are different. The density matrices corresponding 

(k) 

to the right-hand eigenstates are composed from the 
dark eigenstates ) of the Hamiltonian 0, as 



4 4) = ±{\*%W$\ + \*$)(* 



(8b) 
(8c) 



(8d) 



As for the density matrix representation of the left-hand 
eigenstates , the first three g£^ are given by Eqs. I|8a|) 
- ||Sc|, and tho fourth one is g^ = -j^I- The final 
state of the system, after the relaxation has finished, is 



my R 



since in this state 



T^out = 0. By using the density operator representation 

.(*) 



given by r out 

dt ' 

Eq. JSJ of the eigenvectors Tj*L of M, the input-output 
transformation can be written in a simple form as g out = 
Taidin), where T a (g) reads 

T a {Q) = ^ + i(l-Tr m)P D , V = P D gP D (9) 

where Pr> is a projector into the dark subspace of the 

Hamiltonian ©, -Pd = Efc=i I^X^I- 

(/?) XTie case R p ,j e xt > 0: Now the excited state is 
repumped from all external decay channels incoherently 
with a rate of R p . The linear differential equation that 
governs the time evolution of the density operator takes 
the form 4r(r — r) = M'{r — r), where the constant 
vector r satisfies M'r = —d, and the density matrix 
corresponding to r is 

g = sin 2 <p\g+)(g + \ + cos 2 <p\g_) (g- \ 



^ + -»-) sin2( p\g + )( g _ 



| + h.c). (10) 



The left- and right-hand zero-subspaces of the matrix M ' 
coincide and they are three-dimensional. The eigenvec- 
tors rW, [i — 1, 2, 3) satisfy the equation M'(r^ —r) = 
0, and the corresponding density matrices are given by 
Eqs. i|8a[l - l(5c) l. Instead of the mapping in the case (a), 
the input-output states are connected through the rela- 
tion r out = r+^2 i=1 (r^\rii l —r)r^ , which in the density 
matrix representation reads g ut — %{din), where % (g) 
is defined as 

Tb(& = Q-Q +g' + hl-Tr{$})P D , (11) 



where the prime denotes projection into the dark sub- 
space as in Eq. 

Now we turn our attention to finding a pulse sequence 
that yields a desired final density operator of the form 
Eq. 10. The transformation of an initial density oper- 
ator is described by the subsequent applications of the 
mappings of Eqs. © or ljTT|) 



(12) 



where T(g) is equal to T a (g) or %{g). We have to choose 
the number of steps AT, then to find the relative pulse 
amplitudes and phases, defined in Eq. I|T]). by means of 
minimizing numerically the functional 



J({£},Qin,Qf) = (1 -Tr{g f £f} 



1/2 



(13) 



which is the mismatch between the obtained ~Q t (Eq. (JUJ) 
and the required gi (Eq. ©) final density operators. The 
numerical optimization can be performed by means of 
e.g. the conjugate gradient method j^. Due to the spe- 
cial linear property of the mappings T{pipi + P2P2) = 
PiT(pi) + P2T{p2) for pi + P2 = 1, it is sufficient to 
study the convergence for pure initial states, which are 
the arbitrary linear superpositions of the ground states. 
Our aim is to reach a prescribed destination density op- 
erator by applying the same fixed laser pulse sequence for 
all initial states. 

Let us consider a concrete example to demonstrate the 
efficiency of the proposed state engineering method: We 
choose the destination density operator as 



Qi 



+ r } \4 2) )(4 2) \ 



(14) 



with two pure states 
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First we discuss the case when 7 cxt = R p = 0: The ini- 
tial set TL is obtained by discrctizing the four-dimensional 
parameter space - two relative phases and two relative 
amplitudes - describing the possible pure initial states. 
Then we have taken a four-step excitation process, i.e. 
N = 4 in Eq. (fT^|) . and used the conjugate gradi- 
ent method to minimize the functional J{{£ }, di n , gi) of 
Eq. (fT3"|l on the subset TL, g m = |'0in)(^/'in| and iV'in) £ TL. 
The outcome of the optimization is a sequence of four po- 
larization angles and relative phases (ip^ , 9^ l \ fxi} , pf+) 
for I — 1, ... ,4, which characterize the pulse sequence 
{£}. This pulse sequence effects for any initial state such 
a final state, for which the mismatch Eq. (|13(l is less than 
w 10~ 5 (limited by machine-precision). The subsequent 
stages of the transformation of the initial set TL are shown 
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FIG. 2: (Color Online) The transformation of the initial- 
state set after the first- (a), second- (b), third- (c), and 
fourth- (d) excitation steps, shown in the Bloch-sphere of the 
two-dimensional dark subspace. The coordinates are defined 
through the relation g — |(1 + xa x + ya y + za z ), where a q 
are the Pauli's spin operators. The projections of the distri- 
butions to the coordinate planes are also shown. 



in Fig. After the first pulse (Fig. |2fc.) the closure of 
the transformed initial set is the surface of the Bloch 
sphere. The second step (Fig. Eb) yields an elongated 
cigar-shape, while the third step an ellipsoid (Fig. 
distribution, finally the fourth step (Fig. [2Ji) contracts 
the distribution to a point-like region in the Bloch sphere 
with radius around w 10~ 5 . Then we solved numerically 
the Master equation inserting the obtained optimal 
pulse sequence, with a constant Q, — 1 and y m = 1. The 
pulse duration for each step should be chosen so that the 
relaxation process into the actual dark subspace termi- 
nates practically. These times can be estimated from the 
eigenvalues of the matrix M: the one with the smallest 
absolute real part limits the speed of the convergence. 
We note that we have found convergence for any other 
prescribed final state as well, however, the required num- 
ber of steps depends on the purity of the target state: the 



purer the state (i.e. pj> <C p^ or p\ ^ ^> p^ in Eq. 



), the larger the number of steps required. 

For non-zero 7 ex t and R p the optimization of the pulse 
sequence can be done as in the previous case. In a four- 
step process, the numeric optimization yielded a pulse 
sequence for which the mismatch (Eq. is less than 

« 1CU 5 . We have found that the shape of the initial 
distribution transforms in the course of the subsequent 
stages of the excitation process in the same manner as 
before. Then we solved numerically the Master equation 
using the obtained optimal pulse sequence, setting the 
Rabi frequency f2, the decay constants 7i n and 7 ex t, and 
the repumping rate R p to unity: we have found that the 



process converges similarly to the previous case. We note 
that the ratios of 7i n ,7cxt, and R p influence the rapidity 
of the convergence. For sufficiently long time steps the 
process always converges. 

In summary we have worked out a scheme to create any 
pure or a wide class of mixed states in a four-state degen- 
erate A system. Our method is based on an excitation- 
relaxation process, that drives the state of the system 
into the dark subspace of the Hamiltonian that governs 
the dynamics without the decay processes. Although our 
method is not adiabatic, it is robust, because the final 
state is insensitive to the fluctuations in the pulse-area 
of the applied laser-field. A particular advantage of the 
method compared to the adiabatic schemes, that here we 
have a larger freedom to choose the field amplitude. 
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